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ABSTRACT 

We perform a statistical study on the distribution and dynamics of voids in the 2dF- 
GRS. Our statistics are tested and calibrated using mock 2dFGRS catalogues. We 
analyse the dynamics around voids in mock and real 2dFGRS surveys. The void-galaxy 
cross-correlation redshift-space distortions show evidence of the continuing growth of 
voids, confirming another prediction of the hierarchical clustering scenario. A non- 
linear outflow model can be used to provide quantitative estimates of the outflow 
velocities around 2dFGRS voids. These results are consistent with maximum outflows 
of llOkm/s, 210km/s and 270km/s for voids of < r^oid >= 7.5, 12.5 and 17.5h"iMpc, 
assuming a galaxy bias of 6 = 1. As an application for future surveys, our study of 
the mock catalogues shows that direct measurements of the expansion of voids can 
be obtained using peculiar velocity data. We find that it would also be possible to 
detect differences in the velocity dispersion of galaxies in the directions parallel and 
perpendicular to the void walls. 
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. 1 INTRODUCTION 

' The large scale structure in the observed galaxy distribution 
I is characterised by clusters, walls, filaments and also nearly 
1 empty regions called voids (Vogeley, Geller & Huchra, 1994 
and references therein) . The existence of voids has been con- 
firmed by several surveys at a variety of wavelengths (Hoyle 
& Vogeley, 2002), such as the Center for Astrophysics Survey 
(CfA, Vogeley et. al, 1991; 1994), the Southern Sky Redshift 
Survey (SSRS, Gaztafiaga & Yokohama, 1993), the Point 
Source Catalogue Redshift Survey (PSCz, Hoyle & Vogeley, 
2002), the Infrared Astronomical Satellite (IRAS, El- Ad et. 
al, 1997), the Las Campanas Redshift Survey (LCRS, MuUer 
et. al, 2000), the 2degree Field Galaxy Redshift Survey (2dF- 
GRS, Hoyle & Vogeley, 2004; Croton et. al, 2004; Patiri et. 
al, 2005) and the Sloan Digital Sky Survey (SDSS, Hoyle et. 
al, 2005; Rojas et. al, 2005). In all these surveys, voids have 
been found to occupy the greatest volumes in the universe. 

Voids are found in numerical simulations as well, either 
identified using the distribution of dark-matter or haloes. 
Large scale cosmological simulations have been used to 
study the growth and evolution of voids (Regos & Geller, 
1991; Dubisnki et.al, 1993 and references therein). Impor- 
tant progress in understanding the evolution of voids, dis- 
tribution of void sizes and void density profiles was made re- 
cently by several authors using cold dark matter numerical 



simulations and mock galaxy catalogues (for instance Ben- 
son et al. 2003; Goldberg & Vogeley 2004; Sheth & van de 
Wiegaert 2004; Patiri et al. 2004; Shandarin et al. 2004; Col- 
berg et al 2005; Padilla, Ceccarelli & Lambas 2005 and refer- 
ences therein). Other studies have focused on the properties 
of halos populating the inner volumes of voids (Antonuccio- 
Degelou et al. 2002; Gottloeberg et al. 2003; Goldberg & 
Vogeley 2004; Patiri et al. 2004; Goldberg et al. 2005). An- 
other widely used application of void statistics is that of 
probing the bias in the galaxy distribution using large red- 
shift surveys (Mathis & White 2002; Arbabi-BidgoU et al. 
2002; Benson et al. 2003; Goldberg & Vogeley 2004). 

Void statistics have also been found to provide invalu- 
able information on higher order clustering (White 1979; Fry 
1986) which can be used to probe models for galaxy cluster- 
ing (Croton et.al, 2004) and void properties such as sizes, 
shapes and frequency of occurrence, and how these proper- 
ties vary with galaxy type. In general, void statistics may 
provide important clues on the galaxy formation processes 
and can be used to place constraints on cosmological models 
(Peebles 2001). 

Sheth and van de Weygaert (2004) have developed a 
model for the distribution of void sizes and evolution, and 
proved that surviving voids in the Universe should be in ex- 
pansion at the present time. Using simulations in redshift- 
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space, Ryden (1995) found that voids show elongated shapes 
in the direction of the line of sight, a confirmation of the ex- 
pansion of voids. Regos & Geller (1991) have obtain outflow 
velocities from simulations of the growth of voids within 
the large scale structure traced by galaxies. More recently, 
Padilla, Ceccarelli & Lambas (2005) have detected the ex- 
pansion of voids, both from the redshift space distortions 
of void-mass and void-galaxy cross correlation functions in 
numerical simulations, and quantitatively by measuring di- 
rectly positive outflows of void shells. 

Based on these previous results, the central idea in this 
work consists on performing a statistical study of dynam- 
ical properties of galaxies around voids in real and mock 
galaxy catalogues, to prove whether current observational 
catalogues can be used to confirm the predictions from nu- 
merical models. We will develop several statistical proce- 
dures which when applied to mock galaxy catalogues will 
answer these questions. Finally, wo will put the hierarchical 
clustering paradigm to the test once more, by performing 
several void statistics using the 2dFGRS data. 

This paper is organised as follows. We introduce the 
observational and synthetic catalogues in section 2. In sec- 
tion 3 wc describe the algorithm developed to identify voids 
and we apply it to the observational and mock catalogues. 
In section 4 we analise the density profiles of mock galax- 
ies around voids in real and redshift space; and in section 5 
we use redshift space distortions, peculiar velocities, and a 
non-linear approximation to determine properties of the pe- 
culiar velocity field around voids, including the amplitude 
of the expansion of voids and the dispersion of galaxies in 
the directions parallel and perpendicular to the void walls. 
Finally, in section 6, we present a summary of our results. 



2 DATA 

2.1 Observational Data: 2dFGRS 

We use the 2dF Galaxy Redshift Survey (2dFGRS) fi- 
nal data release (Colics et. al, 2003) as our observational 
dataset. The 2dFGRS is one of the widest spectroscopic data 
sets available to date, and contains spectroscopic redshifts 
for approximately 230, 000 galaxies. The source catalogue for 
the 2dFGRS is a revised and extended version of the APM 
galaxy catalogue from which a set of target galaxies, char- 
acterised by extinction-corrected magnitudes brighter than 
bj=19.45, was selected for the construction of the 2dFGRS. 
The main survey regions are two declination strips, plus 99 
randomly located fields. Our results are obtained using the 
two strips described below. One strip is located near the 
southern Galactic pole, spanning approximately 80x15 de- 
grees; the other strip is centred in the Northern Galactic pole 
spanning 75x10 degrees along the celestial equator. The sur- 
vey covers a total of 2, 000 square degrees and has a median 
depth of z = 0.11. 

The completeness of the survey varies with position due 
to unobserved fields, observed objects with poor spectra, 
and objects that could not be observed due to either fiber 
collision constraints or broken fibers. The completeness is 
due in the range [0, 1] . To match the angular selection 
function in the construction of random catalogues we use 



the software developed by the 2dFGRS team and distributed 
as part of the data release that, for any coordinates {a, 5), 
gives the expected probability of galaxy being contained in 
the catalogue. 

2.2 Synthetic Data: Mock 2dFGRS 

We use a ACDM numerical simulation populated with GAL- 
FORM semi-analytic galaxies as our starting point for cre- 
ating mock 2dFGRS catalogues. The numerical simulation 
follows the evolution of 250^ particles in a "Concordance" 
Cosmology with parameters in line with the results from 
the Wilkinson Microwave Anisotropy Map (WMAP, Spergel 
et al. 2003) on the temperature fluctuations of the Cosmic 
Microwave Background, and clustering measurements from 
the 2dFGRS and Sloan Digital Sky Survey (Abazajian et 
al. 2004). Specifically, the box side measures 250h~^Mpc, 
the matter density parameter corresponds to fim = 0.3, the 
value of the vacuum density parameter is f2„ = 0.7, the 
Hubble constant, Hq = /il00s~^km Mpc~l, with h = 0.7, 
and the primordial power-spectrum slope is ns = 0.97. 
The present-day amplitude of fluctuations in spheres of 
8h~^Mpc is set to as = 0.8. The resulting mass per dark- 
matter particle is Mjmn = 1-05 x lO^^h^^Ms. The semi- 
analytic galaxies populating this numerical simulation were 
kindly provided by the Durham group (private communi- 
cation). We briefly explain the procedure followed to this 
end: dark-matter haloes of at least 10 dark-matter par- 
ticles are first identified in the numerical simulation us- 
ing a Friends-of-Friends algorithm with a linking length 
6 = 0.2. The resulting haloes are characterized by a min- 
imum mass of Mmin = 1.05 X lO^^h~'^M0, a median mass 
of Mmin = 2 X lO^^h~^M0, and a maximum mass of 
Mmin — 2 X 15"'^"'^h~^MQ . These haloes are assigned galax- 
ies using the GALFORM code (Cole et al. 2000, Baugh et 
al. 2005). We use these semi-analytic galaxies to construct 
mock 2dFGRS surveys by first positioning an observer in a 
random position within the numerical simulation box, and 
then reproducing the selection function and angular mask of 
the 2dFGRS from this position. This results in a catalogue 
of galaxies of similar properties and observational biases to 
those of the real catalogue. This mock catalogue will be used 
in this work to calibrate our statistical methods, to interpret 
the data, and to detect any systematic biases in our proce- 
dures. In order to do this, we will treat the mock catalogue 
in exactly the same way as the real data. We start our anal- 
yses by defining the same galaxy subsamples for real and 
mock data in the next subsection. 



2.3 Galaxy Samples 

Since the 2dFGRS is an apparent magnitude limited survey, 
the faintest galaxies are only registered at small distances 
while the brightest galaxies cover the whole catalogue. This 
results in a luminosity bias which can be overcome by con- 
structing volume limited samples. 

For a given maximum redshift zum the catalogue is com- 
plete for galaxies brighter than 

Biim = biim - 25 - 5 * log{c * znm) + 5 * log{H), (1) 
where bum is the flux completeness limit and H is the Hubble 
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Table 1. Definition of galaxy subsamples in the NGP and SGP 
regions of mock and real 2dFGRS catalogues. 



constant. This magnitude limit varies from point to point 
in the catalogue, with a brightest magnitude limit cut of 
b < 18.9 (the whole catalogue is complete above this cut). 

We construct several volume limited samples with 
galaxies brighter than Bum and redshifts z < znm- The aim 
to use several samples is is to compare results, differences in 
ziim are related with differences in size of the identified voids 
and differences in Bum are useful to check the underdensity 
stability when using different average galaxy densities. We 
choose 3 limits in redshift and 2 luminosity limits. As a re- 
sult we obtain 12 samples in total for both, the 2dFGRS and 
mock catalogues. 

The parameters and Bii,n used to restrict the sam- 
ples are shown in the second and third columns of Table 
1. The remaining columns of the table show the resulting 
number of galaxies in each subsample for the NGP and SGP 
regions and real and mock catalogues. The values of Bum 
chosen for each subsample satisfy equation 1 for the quoted 
values of zum and bum = 18.9 for 5121, 5151 and 5171, 
bii^ = 19.6 for 5122, 5152 and 5172. 

As can be seen in table 1 the number and number den- 
sity of galaxies diminish when the sample is restricted to the 
brighter galaxies and to higher redshifts. This will have an 
impact on the identification of voids. Different samples will 
be the most appropriate to search for voids of different sizes 
according to the sample volume. For instance, greater vol- 
umes are better suited to find a statistically relevant rmmber 
of large voids. However, in such cases the density of galaxies 
is low and it becomes difficult to find small voids. In the 
opposite case, a smaller volume with a denser population of 
galaxies allows us to identify small voids, but only includes 
a handful of large voids. We take this into account when 
choosing a subsample of galaxies to identify voids of a given 
radius. 



3 VOIDS IN GALAXY CATALOGUES 
3.1 Void finding algorithm 

We construct an algorithm to identify voids in galaxy cat- 
alogues. In order to develop this algorithm we assume that 
voids are approximately spherical. Our assumption is based 
on the approximate shapes of empty regions in numerical 
simulations (Padilla, Ceccarelli & Lambas 2005). The 'a pri- 
ori' knowledge of the void shape suggests the way in that 
the algorithm will work, that is by searching empty spheri- 
cal regions. We define a void as the largest spherical volume 
within which the matter density is below a critical value. Ac- 



cording to this definition, two concentric spheres of different 
radii can not be different voids. Similarly, two or more su- 
perimposed empty spheres are not considered as individual 
voids. 

The algorithm consists of the following steps: First, we 
set a large number of random positions for void centre can- 
didates. These random positions are distributed through- 
out the catalogue. Notice that the efficiency of the algo- 
rithm is partly determined by the density of random centres. 
Therefore, we choose the rmmber of candidate centres such 
that the mean distance between centres is smaller than the 
mean distance between galaxies. Moreover, we always use 
the largest possible number of centres within our computa- 
tional limits. We also check that the relation Vsurvey /Vsphere 
remains large enough for all the void sizes considered. In 
all random positions we consider spheres with radii ranging 
from Tmin to Tmax (cvcry raudom position is the centre of 
a sphere whose radius increases from r,„i„ to Vmax)- Given 
that by definition Tmin is the minimum possible void radius, 
we carefully choose an appropriate value for this parameter. 
The density contrast {5gx) is measured for every sphere. At 
this point the algorithm computes the density within each 
sphere as a function of radius and this value is compared to 
the maximum allowed value, Smax (this allows several candi- 
dates per centre, all with different radii). All the spheres sat- 
isfying the condition 5gx < Smax are considered underdense 
spheres and are selected as void candidates. At this point it 
is still possible that many void candidates share the same 
centre while having different radii. In such a case, we always 
keep the largest sphere. Finally, once the underdense spheres 
have been selected, we decide which underdense spheres axe 
voids. In order to do this we consider three cases: If there is 
a sphere without any superpositions, this sphere is a void. 
If there is a sphere within a larger sphere, the smaller can- 
didate is removed. If there are two or more superimposed 
voids, all of them are removed and the only one considered 
as a void is the largest underdense sphere, which should con- 
tain the others. 



3.2 Voids in mock and 2dFGRS galaxy 
Subsamples 

We apply the void finding algorithm to the mock and 
real 2dFGRS samples defined in the previous section. In 
this work we only identify voids using galaxy positions in 
redshift-space. Padilla, Ceccarelli & Lambas (2005) showed 
that the resulting statistics change only slightly when us- 
ing voids identified in redshift-space with respect to using 
voids identified in real-space. When calculating the density 
contrast required in our algorithm we take into account the 
survey incompleteness. For this aim we use the mask soft- 
ware developed by the 2dFGRS team. 

We find voids with radii ranging from r,nin = 5 Mpc h~^ 
to rmax = 20 Mpc h^^ corresponding to the range allowed 
by the void size constraints suggested by the void probability 
function in the 2dFGRS (Croton et al. 2004; Hoyle & Vo- 
gclcy 2004). Also, according to results obtained by Hoyle & 
Vogcley (2004) void-radii are usually smaller than 25 Mpc/h. 
Another important reason behind the chosen value of rmin 
is given by the galaxy density in our 2dFGRS sub-samples. 
As we are interested in extremely underdense regions we set 
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54 
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2 
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z 
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Table 2. Number of voids identified in the mock and real 2dF- 
GRS, in the NGP and SOP regions. 

Smax = —0.9. This value is in agreement with the mean den- 
sity contrast of voids identified by Hoyle & Vogeley (2004) . 

In table 2 we show the total number of voids obtained 
in each sample. Note that the brightest samples (S083, S103, 
S123, S153 and S173) contain a small number of voids with 
respect to samples constructed using the same z^m and 
fainter galaxies. 

Taking into account the different sizes of voids identified 
in each sample, a more detailed analysis of our samples of 
voids reveals that: 

• In the smallest samples (S121, S122) we do not find a 
significant number of voids with radii larger than lOh^^Mpc. 
However, specially in the most numerous sample (S121), we 
find a good number of voids with radii smaller than about 
7h-iMpc. 

• Samples with the smallest volumes (S121 and S122) are 
more suitable for searching voids of radii between 5 and 10 
Mpc/h rather than larger voids. 

• The intermediate volume samples (S151, S152) give us 
the maximum number of voids with radii of 12, 13, 14 and 
15h~^Mpc. These samples are better suited for searching 
voids of radii between 10 and 15 Mpc/h rather than either 
smaller or larger voids. 

• In the larger volume samples (S171 and S172) we find 
the largest number of large voids {r„oids > 15Mpc/h) 
whereas we do not find any small voids. 

According to this analysis the use of small volume sam- 
ples is suitable for a detailed study of small voids. The final 
selection of samples includes, 
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S121, with z<0.12 and B<-18.88, 5 < rvoid/Mpc/h < 
S151, with z<0.15 and B<-19.36, 10 < Vvoid/Mpc/h < 
S171, with z<0.17 and B<-18.94, 15 < Vyoid/Mpc/h < 



In figure 1 we show the redshift distributions of void 
centres (upper panels) and distributions of void radii (lower 

panels) for voids in samples C121 (solid lines), C151 (dotted 
lines), and C171 (dashed lines), which correspond to z < 
0.12, z < 0.15 and z < 0.17 respectively, for the 2dFGRS 
mock (left panels) and real (right panels) catalogues. As it 
can be seen in figure 1, the distributions are similar in both 
mock and 2dFGRS catalogues, indicating that the general 
characteristics of voids are comparable and therefore the 
mock catalogue is suited to make appropriate comparisons 
to the real data. 
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Figure 1. Radial velocity (upper panels) and void radii (lower 
panels) histograms for voids identified in the (right panels) and 
2dFGRS catalogues (left panels). Different lines correspond to dif- 
ferent samples, solid line for sample S121 , dotted line for sample 
S151 and dashed line for sample S171. 



4 SPATIAL DISTRIBUTION OF GALAXIES 
AROUND VOIDS 

4.1 Galaxy density profiles 

In this section we study the density profiles of galaxies 
around voids. We arc specially interested in the analysis of 
distortions in the profiles caused by the use of galaxy posi- 
tions in redshift-space. We also test the quality of the voids 
identified in this work by examining the density profiles as 
a function of distance to the void centre in the mock cata- 
logues. In all figures, errors arc calculated using the 10 jack- 
knife subsamples of voids (see Padilla, Ceccarelli & Lambas 
2005, for more details on error measurements). In order to 
apply this method for the estimate of errors we follow the 
usual procedure, whereby we calculate the density profile 
for 10 sub- catalogues of galaxies extracted from the sample 
of galaxies under analysis. Each one of these sub-catalogues 
consists of the actual sample with one-tenth of its volume 
removed; this procedure excludes approximately one-tenth 
of the galaxies and each member is only removed once. The 
variance of the results of the true sample is that of the results 
derived from the 10 sub-catalogues. 

The analysis of void profiles resulting from different 
samples generally shows extremely empty voids, with ap- 
proximately zero density of galaxies near the void centres. 
Densities smaller than 5% the median are common in the 
void interiors and the galaxy density abruptly increases at 
distances of about the void radius. This is a good indication 
that we accurately find the boundaries of voids. 

In figure 2 we show the density profiles for voids in the 
mock catalogue in real and redshift space (left panels) and 
for the 2dFGRS survey, only in redshift space (right pan- 
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Figure 2. Galaxy density profiles of mock (left) and 2dFGRS 
(right) galaxies as a function of distance to the void centre. For 
the mock catalogue, we show results in real- and redshift-space 
in gray (cyan) and black lines, respectively. Different line types 
are associated with different subsamples and indicate different 
void radii: solid lines correspond to 5 < r^oid/Mpc/h < 10 and 
sample S121, dotted lines correspond to 10 < r^aid/^V^/h < 15 
and S151, and dashed lines correspond to 15 < r^oi^/Mpc/h < 20 
and SI 71. The sample used for the void identification is the same 
used to trace the profile. 
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Figure 3. Galaxy density profiles as a function of normalized dis- 
tance to the void center. Upper panels: differential density pro- 
files. Lower panels: cumulative galaxy density profiles. Galaxy 
samples and line types are the same as in figure 2 

els). Each figure corresponds to a different sample in the 
mock and 2dFGRS catalogues as is indicated. We analise 
the density profiles as a function of the distance to the void 
centre. As can be seen in this figure the profiles show the 
same behavior regardless of the void radius. However, if we 
anahse in detail the density profiles as a function of void 
radius we can detect some differences. The shape of profiles 
of the smaller voids reach the higher values at a distance of 



Figure 4. Density profiles of mock (left) and real (right) galaxies 
as a function of distance to the void centre. The density profile is 
traced by fainter galaxies than those used to identify the voids; for 
5 < Ti,oid/Mpc/h < 10 (solid line) voids are from S122 and pro- 
files from S121, for 10 < r^oi^/Mpc/h < 15 (dotted line) voids are 
from S152 and profiles from S151, and for 15 < ryoid/Mpc/h < 20 
(dashed line) voids are from S172 and profiles from S171. 



~ 1 void radius. The largest voids only reach profile values 
of ~ 1 at the same distance of 1 void radius. This behav- 
ior is consistent with the results from semi-analytic galaxies 
in the simulation box (Padilla, Ceccarelli & Lambas, 2005). 
Note that the samples used in this figure correspond to dif- 
ferent absolute magnitude limits for different void sizes as 
explained in section 2.3. 

In the case of the mock catalogue, we consider the den- 
sity profiles in real- and redshift-space {Srs and S^s respec- 
tively) and show them in the left panel of figure 2. As can 
be seen in the figure, the profiles show the same shape in 
real and redshift space, although we can appreciate a differ- 
ence at a distance of the order of the void radius where the 
redshift distortion in the void shape introduces an overes- 
timation of the void radius in redshift space. This redshift 
distortion corresponds to the expansion of the void shells, 
which look elongated in the direction of the line of sight 
(Ryden, 1995). This same effect may be responsible for the 
differences in redshift- and real-space profiles at distances 
smaller than the void radius. The density profile in real and 
redshift space tends to be the same at distances greater than 
the void radius. In any case, from the comparison of the rel- 
ative sizes of redshift-space distortion effects and the size of 
error-bars in the density profiles we find that this systematic 
effect is small and therefore we do not correct for it. 

In order to compare with previous works we calculate 
the cumulative density profiles as a function of normalised 
distance to the void centre, and show the results in the lower 
panels of figure 3. The corresponding differential density pro- 
files are shown for comparison in the lower panels of figure 
3. As it can be seen in the lower panels, our cumulative den- 
sities converge to the average density at distances of about 
2 void radius whereas the void profiles by Hoyle & Vogeley 
(2004) remain lower than the mean density. The cumulative 
density slope at the void shell in Hoyle & Vogeley (2004) is 
lower than ours in spite of being selected to have the same 
galaxy density inside the voids. This difference could be a 
consequence of the void identification algorithm. We also 
notice that the profiles shown in the upper panels of figure 
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3, are similar to the results by Benson et al (2003), which 
were obtained from dark matter and galaxies in semianalytic 
simulations. 

We now focuse on fainter galaxies inside voids. In or- 
der to do this, wo calculate the void density profiles using 
fainter galaxies than those used to identify the voids. Figure 
4 shows the results. This analysis is performed in order to 
investigate the possible existence of a faint galaxy popula- 
tion inside voids and also to test variations of the galaxy 
population in voids, and its consequences on void identifi- 
cation and dynamics. As it can be seen in figure 4 there 
are almost no faint galaxies filling the void volumes. We can 
therefore assume that the voids arc indeed virtually empty, 
and therefore the density profile obtained using a complete 
sample of galaxies can approximate the mass density pro- 
file (consistent with semi-analytic simulation results from 
Padilla, Ceccarelli & Lambas 2005). This assumption has 
important dynamical consequences as will be seen below. 



5 DYNAMICAL PROPERTIES OF GALAXIES 
IN VOIDS 

We analyse the dynamics of galaxies around voids following 
two different approaches. The first is a qualitative approach, 
where we study the shape of the isocontours of the void- 
galaxy cross-correlation function measured in the directions 
parallel and perpendicular to the line of sight. From this 
study we can infer characteristics of the peculiar velocity 
field. The second approach is a quantitative measurement 
of galaxy outflows or void expansion, which can be derived 
in two ways; i) the first measurement of the void expansion 
can be achieved by applying a simple, non-linear theory ap- 
proximation that relates the peculiar velocity outfiows to 
the void-galaxy cross-correlation function; this method can 
be applied to the mock and real 2dFGRS data, ii) The sec- 
ond method is only possible in the mock catalogue due to 
the current observational limitations, where we directly cal- 
culate the outflow amplitudes from peculiar velocities. The 
use of peculiar velocities allows a direct, quantitative mea^ 
surement of the outflow of galaxies from voids, which can be 
compared to predictions from non-linear theory approxima- 
tion and results from the rmmerical simulation. We also test 
whether a mock peculiar velocity survey can detect subtle 
properties of the dynamics of galaxies around voids, such as 
a larger velocity dispersion in the direction of the void walls. 

5.1 Dynamics using redshift space distortions 

We have performed a study of the void-galaxy cross correla- 
tion function, that is the correlation function between void 
centres and galaxies, as a function of the directions parallel 
and perpendicular to line of sight, tt and a respectively. 

The assumption that our Universe is isotropic leads 
us to expect that the correlation function should also be 
isotropic. However, if wo measure the correlation function 
in redshift-space instead of in real-space, anisotropies ap- 
pear. This deviations from isotropy are known as redshift- 
spaco distortions. Those distortions are due to peculiar ve- 
locities, and therefore can provide information on the dy- 
namical properties of galaxies around voids, ^(cr, tt) is calcu- 
lated using the estimator 



C(cr,7r) = DD(Aa, A7r)/i3i?(Aa, Att) - 1 (2) 

where DD is the number of galaxy and void-centre pairs in 
the range of distances (Act, Att), and DR is the number of 
void-random galaxy pairs in the same range of distances 

In a hierarchical scenario of structure formation, the ob- 
jects in underdcnse regions tend to be accreted toward more 
massive systems. According to this, the galaxies inside voids 
and in void shells, would be attracted to the matter outside 
of the void, producing a large scale streaming motion in the 
form of an outflow from the void centre or a void expansion. 
The outflow motion would be seen in the distortion pattern 
of the z-space correlation function as an elongation in the 
direction of the line of sight (tt direction). The upper pan- 
els of figure 5 show the contour lines for the samples S121 
from the mock (left panels) and 2dFGRS (right panels) cat- 
alogues, for small void radii (5 < r-^oid/ Mpc/h < 10). The 
two panels in the middle of figure 5 correspond to intermedi- 
ate void sizes 10 < rvoid/Mpc/h < 15 obtained firom sample 
S151. The bottom of figure 5 shows the distortions for the 
largest voids (15 < r^oid/Mpc/h < 20) obtained using sam- 
ple S171. As is mentioned in section 3, the different samples 
considered in the upper, middle and lower panels of the fig- 
ure have been chosen to have the most convenient number 
of voids for a high statistical significance and to reduce the 
noise in the results. The estimated jackknife uncertainties in 
5((T, tt) arc ~ 10% of the value of ^. There is a clear elongated 
distortion pattern in all panels of figure 5 which is consis- 
tent with an outflow motion. The fact that the 2dFGRS 
galaxies show such an elongation around voids represents a 
further confirmation of the hierarchical clustering scenario 
of structure formation. However, the noisy results restrict 
the quantitative applications of this method. Therefore, a 
different approach should be followed in order to quantify 
the expansion of voids in the Universe. 

5.2 Quantitative measures of galELxy outflows 
Eiround voids 

In order to study the statistical characteristics of the pecu- 
liar velocity field around voids we estimate the mean ampli- 
tude of the velocity field as a function of the distance to the 

void centre. We place our reference system at the void cen- 
tres so that a negative velocity indicates an outfiow motion 
(and a positive velocity indicates an infall onto the void cen- 
tre). We adopt two approaches, an indirect derivation of the 
velocity field from a theoretical model, and a direct mea- 
sure of the peculiar velocity amplitude. The first method, 
requires many assumptions about the adopted cosmological 
model in order to derive the velocity field. However, this 
method works satisfactorily even when assuming a simple, 
non-linear approximation. The second method consists of 
the calculation of the velocity field amplitude using galaxy 
pecuhar velocities. Even though this is the most simple and 
natural approach, it has the disadvantage of being only ap- 
plicable to mock catalogues. 

5.2.1 Non-Linear Model for the Outflow 

We apply a non-linear theory model by Croft, Dalton & Ef- 
stathiou (1999) to describe the dynamical behavior of par- 
ticles around underdense regions. We choose to use a non- 
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Figure 5. Redshift space void-galaxy correlation function ^(tj, tt) estimated from mock (right panels) and real 2dFGRS (left panels) 
catalogues, a and tt represent separations parallel and perpendicular to the line of sight. Each panel corresponds to different ranges of 
void radii and are obtained from different samples; top: 5 < r^ai^jMpclh < 10 and S121, middle: 10 < r^aid/Mpc/h < 15 and S151, 
bottom: 15 < Tyoid/Mpc/h < 20 and S171. Dashed lines represent ^ = —0.6,-0.5,-0.4, dotted lines represent t; = —0.3, —0.2, —0.1, 
solid line represents the iso-{ contour corresponding to ^ = 0., dot-dashed lines represent ^ = 0.1,0.2 and short-dashed-long-dashed lines 
represent ^ = 0.3, 0.4, 0.5, 0.6, 0.7. 
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linear model in favor of a linear approximation since in the 
inner regions of voids, the overdensity is significantly dif- 
ferent from zero and therefore, the linear theory condition 
\5\ « 1 is not satisfied. The outfiow model, which is analo- 
gous to the infall model for overdensities, assumes spherical 
symmetry. When analyzing both, mock catalogues and real 
data, we assume the concordance ACDM cosmology consis- 
tent with results from WMAP (Spergel et al. 2003). As the 
overdensity within regions around voids is ~ —1, we can ex- 
pect the non-linear approximation to describe the velocity 
field accurately, 



1/ 



in"'«i/or5(r)/(l + 5(r))''-^5 



(3) 



where 5{r) is the matter density contrast within a sphere 
of radius r. We apply the outflow model to predict the 
amplitude of peculiar velocities. The model uses the mea- 
sured values void-galaxy cross correlation function in red- 
shift and real space (^(s) and ^{r), respectively) to calcu- 
late the galaxy density contrast (Sgj, (r) (Croft et. al, 1999); 
the subindex gx indicates that we can not measure the void- 
mass cross-correlation function directly. Therefore, the infall 
velocities obtained from this formalism actually correspond 
to 



mf 



'■ bVinf 



(4) 



400 



200 



valid where \5\ « 1 (that is, at the void boundaries or 
further away), b is the unknown bias factor relating the mass 
and galaxy density contrasts around voids, which in this 
work is assumed to be b = 1 (see for instance Verde et al. 
2002). 

We first apply this formalism to the void-galaxy cross 
correlation function obtained from the mock catalogue and 
show the predicted peculiar velocities in gray lines in figure 
6, where each line-type corresponds to different void radii 
(same ranges in black lines). 

We then apply the outfiow model to the 2dFGRS cat- 
alogue. We use the redshift space correlations obtained for 
the observational catalogue assuming that our results will 
not change significantly when using redshift-space correla- 
tions as opposed to real-space data. The predicted velocities 
are shown in black lines in figure 6. For comparison, we also 
show the predicted velocities for the mock catalogue in gray 
lines. The maximum outfiow amplitudes for the 2dFGRS, 
and distances where this maximum occurs, are in excellent 
agreement with the predictions for the mock catalogue. Re- 
gardless of the void size, both curves exhibit very similar 
amplitudes and shapes, specially at distances of the order of 
one void radius. 

It is very reassuring to find the same velocity outfiows 
in mock and observational catalogues. This suggests we are 
obtaining a reliable estimate for the velocity of void shells in 
the 2dFGRS. Only at large separations from the void cen- 
tres the outflow velocities became considerably larger in the 
2dFGRS. It remains to be confirmed whether this result is 
to be expected from the data or simply produced by sys- 
tematic effects perhaps induced by the use of redshift-space 
correlations. 
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Figure 6. Predicted velocity field amplitude as a function of 
distance to the void centre for the 2dFGRS (black lines) using 
the non-linear velocity model. Different line-types correspond to 
different void radii, and the subsamples shown are the same as in 
figure 7. Predicted velocities for the mock catalogue are shown in 
gray lines for the same samples. 



5.2.2 Galaxy outflows from mock peculiar velocities 

Studies on the outflow of void shells have been performed 
mainly in numerical simulations and semi-analytical galaxies 
by several authors. For instance, Regos & Geller have stud- 
ied the evolution of voids in numerical simulations obtain- 
ing the peculiar streaming velocities of void walls (Regos & 
Geller 1991). Dubinski et al. (1993) have analised the pecu- 
liar velocity field surrounding voids in simulations. Another 
related work was performed by Seth & van de Weygaert 
(2004) who studied the void size evolution in simulations. 

We compute the amplitude of the velocity field around 
voids using the line-of-sight component of galaxy peculiar 
velocities in mock 2dFGRS catalogues, assuming a radial, 
spherically symmetric outfiow motion. Note that as the real 
2dFGRS survey contains no data regarding peculiar veloci- 
ties, this method can not be applied to this particular cat- 
alogue. Therefore, all the analyses in this subsection are re- 
stricted to the mock catalogues. 

In Figure 7 we show the peculiar velocity outfiow as 
a function of the distance to the void centre, derived from 
the mock peculiar velocities (black lines, the model results 
are shown in gray lines for comparison). Different line types 
correspond to different void radii and samples. As can be 
seen in the figure, the minimum velocity occurs at distances 
of about the void radius, so that larger voids show a mini- 
mum velocity amplitude (maximum outflow) at greater dis- 
tances from the void centre. The absolute velocity minimum 
is larger for greater voids and the range of maximum outfiow 
velocities goes from about 150fcm/s for the smallest voids 
to more than 200fcm/s for the largest voids. By comparing 
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Figure 7. Outflow as a function of distance to the void cen- 
tre derived from peculiar velocities in tlie mock catalogue (black 
lines). Different lines correspond to different ranges of void ra- 
dius and are obtained from different subsamples (solid lines, 5 < 
r^aid/Mpc/h < 10 and S121 , dotted lines, 10 < Tyaid/^'^pc/h < 
15 and S151, dashed lines, 15 < r„M/Mpc/h < 20 and S171). 
Gray lines show the velocity amplitude predicted by the outflow 
model using the same samples and void radii as for the black lines 



Figure 8. Velocity field amplitude as a function of the distance 
to the void centre when including distance measurements errors. 
Dark-gray lines without error bars correspond to the measured 
outflow velocities using peculiar velocity data with no errors; 
light-gray lines show the results when including errors in the pecu- 
liar velocity; black lines show the results when only using peculiar 
velocities (with errors) of galaxies within 45 degrees of the line 
of sight from the void centres. Different line types correspond to 
different ranges of void radius and samples. 



between measured and model outflow velocities (black and 
gray lines, respectively) we can notice a qualitative agree- 
ment, specially for the larger voids. Small voids show some 
discrepancies at distances a few times the void radius; this 
can be due to failure of the approximations assumed in the 
identification of voids of small size, such as spherical sym- 
metry. 

The maximum measured outflow velocities are success- 
fully described by the model, as well as the distance to the 
void centre where the maximum outflow velocity occurs. 
Inspection of figures 6 and 7 shows that model and mea- 
sured velocities around medium and large size voids (dot- 
ted lines and dashed lines, respectively) first reach a min- 
imum value and then increase in value to eventually reach 
Vfiow ~ Okm/s. Notice that for the largest voids there are 
still negative velocities at large distances. Small voids (solid 
lines) show positive velocities at distances of about two void 
radii, and stay on positive values over a wide range of dis- 
tances to finally return to Vfiow ~ km/s. This behavior 
could be related with an infall of void shells due to surround- 
ing overdensities. It is also related with the evolution of voids 
analised by Sheth and van de Weygaert (2004) in large scale 
numerical simulations, where they use two processes to de- 
scribe the void evolution: void-in-void processes and void- 
in-cloud processes. The void-in-cloud refers to small voids 
embedded in dense regions, which show collapse instead of 
expansion. Then, the infall of small void's shells can be re- 
lated to the collapse in the void-in-cloud process. 

In order to test the effect of distance measurement er- 



rors on the velocity field amplitude, we apply a random un- 
certainty proportional to 10% of the redshift-space distance. 
As expected, the main effect of distance uncertainties is to 
overestimate the peculiar velocity field amplitude specially 
in the deepest samples. The resulting velocity amplitudes 
are shown in light-gray lines in figure 8 whereas in dark-gray 
lines we show the resulting velocity outfiows when distance 
measurement errors are not included. 

As can be seen in figure 8 there is a strong shift of out- 
fiow velocities to more negative values. For example for the 
greatest voids (dashed lines), the maximum outfiow velocity 
changes from 200 km/s to 400 km/s, a systematic increase 
of 100% when including distance measurement errors. 

The largest contribution from distance measurement er- 
rors to the systematic offset in the outfiow velocity comes 
from using radial peculiar velocities to infer the outflow ve- 
locity when the angle subtended by the observer and the 
galaxy from the void centre is close to 7r/2. In this case, the 
component of the peculiar velocity in the radial direction 
to the void centre is nearly null, and the deprojection pro- 
duces large spurious outflow velocities responsible for the 
systematic offset. Therefore, we repeat the analysis using 
only galaxies within 7r/4 of the line-of-sight (as seen from 
the void centres) and show the results in black lines in fig- 
ure 8. As can be seen the systematic offset is much less 
important in this case. 
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5.2.3 Maximum outflow amplitude 

The results obtained in the previous subsections strongly 
suggest that there is a relation between the the void size and 
the distance where the maximum outflow occurs. The void 
radius seems to determine in a very simple way the distance 
of the maximum outflow velocity. In order to characterise 
in a more quantitative way such a relation we study the 
contour lines of outflow velocity as a function of void radius 
and distance to the void centre. 

The contour lines of outflow velocity directly measured 
from the mock peculiar velocities are plotted in the left panel 
of Figure 9, where the x-axis represents the void radius and 
the y-axis the distance to the void centre. The solid lines 
in figure 9, which represent positive velocities (i.e. galaxy 
infall to the void centre), are only found for small voids 
(r-uoid < llMpc/h) at relatively large distances compared 
to the void size. It is easily seen how large outflow veloci- 
ties {vout > lOOkm/s) are found around the largest voids 
(with radii r^oid > 14:Mpc/h) at distances in the range 
10 < rMpc/h < 20. Note that this range is approximately 
centred at the average void radius. The right panel shows 
the results when adding distance measurement errors to the 
mock peculiar velocities. As can be seen, there is also a clear 
relation between void radius and the distance corresponding 
to the maximum outflow velocity; however, as a consequence 
of large velocity errors, the resulting outflow velocities are 
much higher in this case. 

In order to characterise the relation between the void 
radius an the distance at which the minimum velocity oc- 
curs, we flnd that a linear relation fits the results: 

r{vmin) = air-yoid + bi, (5) 

where ai = 0.7 and bt — 2. This expression is represented by 
a gray dotted line in both panels of figure 9. Although this 
is the simplest way to describe the relation r{vmin)-rvoid we 
also fit an exponential relation: 

r-(umin) = expir^oidjaf (6) 

which is shown in a gray solid line in figure 9, with best 
fitting parameters a — 0.62 and b — 0.30 (no errors) and 
a = 0.55 and b = 0.33 for the results when including dis- 
tance measurement errors. This exponential expression can 
also fit the distance where the velocities reach 0km/ s; in this 
case, the best fitting parameters are a = 1 and b = 0.5 (when 
not considering errors) . The parameter values obtained from 
peculiar velocities with no errors are similar to those ob- 
tained by fitting the same expression to semianalytic galax- 
ies (Padilla et al, 2005). As can be seen in figure 9, the linear 
and exponential functions are very similar, and both can be 
used to give an approximate relation between radius and 
maximum outflow amplitude distance. 

We analise the dependence of the maximum outflow dis- 
tance on void radius predicted by the non-linear outflow 
model, in this case for both mock and observational cata- 
logues. The aim in doing this is to compare the results from 
both catalogues, and to test the quality and efficiency of 
the model to describe the velocity fleld. We produce the ve- 
locity contour lines predicted by the model and show them 
in figure 10, where the left panel corresponds to the mock 
catalogue and the right panel corresponds to the 2dFGRS. 

If we compare both panels of figure 10 we note that 




Figure 9. Contour lines of outflow velocities as a function of 
void radius (x-axis) and distance to the void centre (y-axis) 
for the mock catalogue, using peculiar velocities. The left panel 
shows the results from peculiar velocities with no errors; the 
right panel shows the results when distance measurement er- 
rors are added. Dot-dashed lines represent velocities lower than 
— 160km/s (—200, —180, —160), long-dashed lines represent veloc- 
ities -160 < v/km/s < -100 (-140, -120, -100), short-dashed 
lines represent velocities —100 < v/km./s < —40 (-80, —60, —40), 
dotted lines represent velocities —40 < v/km/s < (—20, —10, 
—6, —2, —1), and solid lines represent positive velocities (0, 5, 15, 
20). 

from a qualitative point the results from mock and theory 
are quite similar. In a more quantitative analysis, we find 
that the maximum velocities around large voids {r^oid > 
14h-^Mpc) occur at scales in the range 10 < rh"^Mpc< 20 
in both cases. The results shown in this figure are also in 
agreement with those found from the direct measurement of 
the peculiar velocity field shown in figure 9. 

We fit the maximum model outflow velocities using 
equations 5 and 6. The results from the mock catalogue are 
best fitted using the exponential expression with a = 0.90 
and b = 0.33; the results from the 2dFGRS theory predic- 
tions are fitted by a = 1.05 and b — 0.35, also for the ex- 
ponential fit. Note the similarity between the results from 
mock and real data. 



5.2.4 Velocity dispersion in void shells 

In order to fully characterise the dynamics of galaxies in void 
shells, we also compute the peculiar velocity dispersion in 
the directions parallel and perpendicular to the void shell. 
Again, due to the lack of peculiar velocity information in 
the 2dFGRS survey, all the analyses in this subsection are 
restricted to the mock catalogues. 

As we only have the line-of-sight projection of the pe- 
culiar velocity, the best signal for velocities perpendicular to 
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Figure 10. Contour lines of peculiar velocity outflows derived 
from the theoretical model, as a function of void radius (x-axis) 
and distance to the void centre (y-axis), for the mock and 2dF- 
GRS catalogues (left and right panels respectively). Gray solid 
lines represent velocities greater than —340 km/s (—340, —320, 
-300, -280), gray dotted lines -280 < v/km/s < -200 (-260, 
-240, -210), black dot-dashed lines -220 < v/km/s < -140 
(-200, -180, -160), black long-dashed lines -160 < v/km/s < 
-100 (-140, -120, -100), black short-dashed lines -100 < 
v/km/s < -40 (-80, -60, -40), black dotted lines -40 < 
v/km/s < (—20, —10, —6, —2, —1), black solid lines represents 
positive velocities larger than 10 km/s (0, 5, 15, 20), cyan short- 
dashed lines 20 < v/km/s < 100 (40, 60, 80), cyan long-dashed 
lines 160 < v/km/s < 100 (100, 120, 140), and cyan dot-dashed 
lines V > 160fcm/s (160, 180, 200). 

void shell (in the void radial direction) comes from galaxies 
close in projection to the void centre, whereas for the veloc- 
ity parallel to the void shell, this comes from galaxies away 
from the void centre, in projection. Taking this into account 
we use galaxies within 45 degrees of the line of sight from the 
void centre when calculating the dispersion in the direction 
perpendicular to the void shell, and only galaxies at more 
than 45 degrees with respect to the line of sight from the 
void centre when calculating the dispersion in the direction 
parallel to the void shell. 

Figure 11 shows the relative velocity dispersion as a 
function of distance to the void centre in the mock catalogue. 
Black lines show the dispersion about the mean motion in 
the radial direction, and gray lines show the dispersion in 
the direction parallel to the void walls. Solid lines correspond 
to small voids from sample SlOl, dotted lines correspond to 
intermediate size voids from sample S151 and dashed lines 
correspond to large voids from sample S171. Errors are cal- 
culated considering both Poisson and jackknife errors, added 
in quadrature. 

It can be seen in figure 11 that the dispersion is not 
very dependent on distance to the void centre, specially for 
the smallest voids. The plot also shows that parallel veloci- 



Figure 11. Galaxy peculiar velocity dispersion relative to the 
void shell as a function of normalised distance to the void cen- 
tre. Black lines show the dispersion in the direction perpendic- 
ular to the void shells and gray (cyan) lines indicate dispersion 
in the direction parallel to the void shells. Different line-types 
correspond to different void radii and galaxy samples: solid lines 
are for 5 < r^^id/Mpc/h < 10 and S121 , dotted lines are for 
10 < r^aid/^pc/h < 15 and S151, and dashed lines are for 
15 < r^oid/Mpc/h < 20 and S171. 



ties are always higher than radial velocities indicating that 
galaxies tend to move faster in the direction of void shell. It 
can also be seen that the smallest voids show larger values of 
velocity dispersion than the largest voids. All these results 
are consistent with those obtained using semianalytic galax- 
ies in a full simulation box (Padilla, Ceccarelli & Lambas, 
2005). 

We also show in figure 12 the ratio between the velocity 
dispersion in the directions parallel and perpendicular to the 
void shell. In order to improve the statistics we have used 
a combined sample of intermediate and larger voids 10 < 
rvoid < 20. For this sample the ratio between parallel and 
perpendicular dispersion is shown in figure 12. In the inner 
panel are shown the ratios (Jparaiiei/o-perpendicuiar obtained 
from the same samples analysed in figure 11 

As it can be seen in the large panel of figure 12, the ratio 
of parallel to perpendicular dispersion is nearly constant at 
(^parallel /(^perpendicular ~ 1-5 for distauces to the void Centre 
in the range 0.5 < r/r^oid < 2.5. Note that these values 
are comparable to those found in numerical simulations by 
Padilla, Ceccarelli & Lambas (2005). 

We add distance measurement errors to the galaxy pe- 
culiar velocities proportional to the distance to the observer 
(see for instance Padilla & Lambas, 1999). The velocity dis- 
persions obtained using peculiar velocities affected by dis- 
tance measurement errors are shown in figure 13, where 
black lines show the velocity dispersion along the direction 
perpendicular to the void shell and gray lines show the dis- 
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Figure 12. CTparallel/<^p erpendicular a function of normalised 
distance to tlie void centre for 10 < r^„i^/Mpc/h < 20; in- 



ner panel: (t„ 



parallel / ^perpendicular 



as a function of normalised 



distance to the void centre. Different line-types correspond 
to different void radii and galaxy sample: solid lines corre- 
spond to 5 < r^„ici/Mpc/h < 10 and S121 , dotted lines to 
10 < r^aid/Mpc/h < 15 and S151 and dashed lines to 15 < 
rvoid/Mpc/h < 20 and S171. 



Figure 13. Galaxy peculiar velocity dispersion relative to the 
void shell as a function of normalised distance to the void centre 
affected by distance measurement errors. Black lines indicate the 
dispersion in the direction perpendicular to the void shell, and 
gray (cyan) lines indicate the dispersion in the direction parallel 
to the void shell. Different line-types correspond to different void 
radii and galaxy samples as in figure 11. 



persion parallel to the void shell. Solid lines correspond to 
small voids from sample SlOl, dotted lines correspond to 
intermediate voids from sample S151 and dashed lines cor- 
respond to large voids from sample S171 as in figure 11. 

Inspection of figure 13 shows that velocity dispersions 
are higher when peculiar velocity errors are included. We can 
also see that this effect is more important for the largest 
voids. The reason behind this systematic offset is that we 
use samples that reach higher redshifts for larger voids. In 
spite of this, a qualitative comparison of figures 11 and 13 
shows that the behavior of the velocity dispersion parallel 
and perpendicular to the void shells are similar in both plots 
and that it can still be detected that galaxies move faster 
along the void walls, or alternatively, filaments. 



6 CONCLUSIONS 

In this work we have studied statistical properties of the 
dynamics of voids and galaxies populating void shells. We 
were specially interested in detecting and characterising the 
expansion of voids, and in order to do so, we studied the 
redshift space distortions detected in the galaxy-void cross 
correlation function measured in the directions parallel and 
perpendicular to the line-of-sight. We also analysed the am- 
plitude of the peculiar velocity outfiows obtained from the 
peculiar velocities in mock catalogues and from non-linear 
theory approximation (Croft, Dalton & Efstathiou, 1999). 
In addition, we analysed the velocity dispersion in the di- 



rections parallel and perpendicular to the void shell. In all 
our statistical analyses we test the effect of peculiar velocity 
errors in our results. We summarise our results as follows: 

• We applied a void finding algorithm, which is capable 
of accurately identifying voids in galaxy catalogues affected 
by angular and radial selection functions, to volume limited 
samples of mock and observational 2dFGRS catalogues. 

• The redshift space distortion of galaxies around voids 
shows elongations along the line of sight. This can be inter- 
preted as evidence of an outfiow motion in both mock and 
observational data. 

• Outflow velocities from the mock catalogue (without 
errors) are consistent with the non-linear theory predictions 
based on measurements of the correlation functions in the 
mock catalogue. 

• By applying the non-linear outflow model to 2dFGRS 
data, we flnd that larger voids would be characterised by 
larger maximum outflow velocities: llOkm/s, 210km/s and 
270km/s for < r^oid >= 7.5, 12.5 and 17.5h-^Mpc respec- 
tively; assuming a galaxy bias b — 1. We flnd that out to a 
distance of the order of the void radius. The results from the 
application of the non-linear outfiow model to observational 
data are consistent with results from the mock catalogues . 

Our analysis of mock catalogues indicates that several 
additional results on the dynamics around voids could be 
derived from the 2dFGRS survey if peculiar velocity data 
was available as well, or at least for a few tens of thousands of 
galaxies sampling a large enough volume so that a reasonable 
number of voids could be found. 
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• By studying the mock catalogues, we are able to per- 
form a direct measurement of galaxy outflows using peculiar 
velocities. Outflows of about 200 km/s are measured, in ac- 
cordance with the results from numerical simulations from 
Padilla et al.(2005). 

• When distance measurements errors are included in the 
mock galaxies, the estimated outflow velocity is biased to 
higher values, which could be up to 400 km/s for the larger 
voids. 

• Measurements of the velocity dispersion in shells 
around voids in the mock catalogues show that these are 
larger in the direction parallel to the void walls, in agreement 
with results from the full numerical simulation (Padilla. Cec- 
carelli & Lambas, 2005). This result indicates that galaxies 
tend to move preferentially in the direction parallel to void 
shells, that is, along large-scale structure filaments and walls, 
and that a peculiar velocity catalogue with sufficient data 
could allow this measurement to be done, constituting a new 
test to the hierarchical clustering paradigm. 
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